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Abstract 

We consider heavy quark interactions in quenched and unquenched lattice QCD. In a region just above the 
deconfinement point, non-Abelian gluon polarization leads to a strong increase in the binding. Comparing 
quark-antiquark and quark-quark interaction, the dependence of the binding on the separation distance r 
is found to be the same for the colorless singlet QQ and the colored anti-triplet QQ state. In a potential 
model description of in-medium J behavior, this enhancement of the binding leads to a survival up to 
temperatures of 1.5 Tc or higher; it could also result in J/ijj flow. 



1 Introduction 

The interaction of a heavy quark-antiquark {QQ) pair in strongly interacting matter has been studied in 
finite temperature lattice QCD in the quenched approximation as well as for the cases of two light and two 
light plus one heavy quark species [1-4]. In all these studies, one obtains the difference F(r, T) between 
the color singlet free energy with and without the heavy quark pair, as function of the temperature T of 
the medium and the separation distance r of Q and Q. Schematically, we can write 

Fo(T) = -T\nj dTexp{-no/T} (1) 

and 

FQir,T) = -T In j exp{-7^Q(r)/T} (2) 

for the free energies in question; here the Hamiltonian Tig describes the interacting quark-gluon plasma 
alone, and Hq^r) that of the plasma containing the static QQ pair; the integral / dT denotes the grand 
canonical phase space integration and summation. The color singlet free energy difference provided by 
lattice studies is then defined as 

F{r,T) = FQir,T)-FoiT). (3) 

In Fig. [U we show the r-dependence of F(r,T) at different temperatures above the deconfinement point, 
i.e., for T > Tc- The vacuum form (for T = 0) also shown here is the Cornell potential [5] 

V{r) = a r --, (4) 
r 

defined in terms of the string tension a and the Coulomb coupling a ~ 7r/12. For very small r <C T~^, the 
small color singlet pair is neither seen nor affected by the medium, so that the interaction is specified by 
the T-independent running coupling a(r), which for larger r saturates at the canonical string value 7r/12. 
On the other hand, at very high temperatures and comparatively large r ^ T~^, in the perturbative 
limit, we expect an r- independent running coupling a(T). 

In full QCD below Tc, the presence of light quarks (g) leads to qq pair production and string breaking; 
from quarkonium studies, the string breaking energy at T = is estimated to be about 1.0 - 1.2 GeV. 



In the quenched case below Tc, the free energy diverges in the large distance limit, with a temperature- 
dependent string tension cr(T) which vanishes at T = Tc. Temperatures above Tc, however, as seen in 
Fig. [T^, also lead to a finite large-distance limit, with a temperature-dependence which is very similar 
to that found in full QCD. Since here there are no light qq pairs to provide string-breaking, the large 

^Based on joint work with O. Kaczmarck and F. Karsch 
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Figure 1: Free energy difference for a color-singlet 
two-flavor QCD [f]. 



pair as function of r, (a) for quenched and (b) for 



distance behavior in the quenched case must arise from gluon screening effects, i.e., it is of a non-Abelian 
origin. We therefore expect that also in full QCD gluonic screening plays a crucial role. 

From -F(r, T) one obtains through standard thermodynamic relations the corresponding difference for the 
internal energy U{r, T), 

Uir, T) ^ -T- (ffi^) = Fir^ T) ^ T (^^) = Fir, T) + T Sir, T), (5) 



and for the entropy S'(r, T), 



5(r,T) 



fdFir, T) 



(6) 



V dT 

The internal energy behavior is shown in Fig. [5] [3,4]. From eqs. ([T|), ([3]) and (O, we obtain 

C/(r,T) = (HQ(r,T))-(Wo(r)), (7) 

which for a static QQ pair, with no kinetic energy, measures the change in potential energy due to the 
introduction of the pair. It is seen in Fig. that just above the deconfinement point, the potential is 
much stronger than at T = 0. In the following, we want to study the origin of this increase. 

Consider first the underlying medium, a plasma of unbound (but interacting) gluons and (for full QCD) 
quarks at temperature T. We had seen that gluon screening plays a decisive role in the medium, and 
when we speak of constituents, this should be kept in mind. Any interaction is mediated by gluons, and 
when the gluonic interaction range is reduced through color screening, this holds for quark interactions as 
well. We denote by n(r) the average density of constituents, so that d = rT^I"^ is the average separation 
distance between adjacent color charges. The interaction range for the constituents is given by the 
correlation function, 

fir) ~ exp{-r/^(r)}, (8) 

which measures the interaction strength between constituents separated by a spatial distance r. The 
correlation length ^(T) is a basic property of the medium, and it will diverge for T — s- Tc in case of a 
second order confinement/deconfinement transition. In the high temperature limit, we expect a non- 
interacting gas and hence ^ — > 0. 

We now add a static QQ pair to our system. As long as the Q and the Q are sufficiently far apart, 
their color charges are effectively screened and asymptotically, they do not interact with each other. 
Nevertheless, the Q as well as the Q individually interact with the medium, and around each charge, 
this will lead to polarization effects in the region where its presence is felt. Hence in a sphere of average 
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Figure 2: Internal energy difference for a color-singlet QQ pair in two-flavor QCD, (a) as function of r 
and (b) for r ^ oo, compared to the corresponding free energy (solid line) [4]. 



radius £,{T) around the Q (and similarly around the Q), there will be an excess of color compared to the 
state of the medium before insertion of the static charges. The polarization screens each static charge, 
so that its strength is reduced the further away we measure it. This evidently is a non-Abelian effect, 
since an electric charge cannot be screened by photons, and in quenched QED the QQ interaction would 
remain purely Coulombic even in the large distance limit. When the QQ separation distance r is reduced, 
eventually interaction between the charges becomes significant. This interaction initially has a Debye- 
screened Coulomb form, with Oefi ~ a exp{— /xr} as the effective coupling, where /i denotes the screening 
mass. It has to be emphasized that such a picture makes sense only for r ^ ^, so that we really have two 
well separated polarization spheres, each of which contains a static charge surrounded by many medium 
charges, with an overall excess of the opposite charge. 

When the heavy quarks are brought still closer to each other, the polarization spheres begin to overlap [4]. 
The overlapping spheres continue to attract each other, but this interaction now has two components 
and can no longer be described by a pure Debye-screened potential. On one hand, there is the direct 
interaction between the Q and the Q, which in the limit of small r becomes just the Coulombic vacuum 
form a/r. On the other hand, the charged constituents inside the overlapping polarization spheres also 
attract each other directly, with a strength determined by the deviation of the region from the state 
without a QQ pair. Once the spheres no longer overlap, this interaction becomes by Gauss' law just the 
Debye-screened Coulomb form with a reduced effective charge. 

In order to determine the behavior of QQ binding in a deconfined medium, we shall first consider the 
large and small distance limits of the thermodynamics potentials. 

2 The Large Distance Limit 

In the limit of large QQ separation, for r — > cx), we have two fully screened and hence non- interacting 
charges. Nevertheless, the thermodynamic potential differences F, U and S do not vanish: they specify 
the effect of the interaction of each of the two independent charges with the medium. This is not related 
to the interaction of Q and Q, as is easily seen by considering the corresponding quantities for a pair of 
heavy quarks (a "diquark"). Lattice calculations [6,7] have in fact shown that for r ^ cxd one obtains 

PQUT)=F§l{T)^2F,{n (9) 

where F'^^^ corrresponds to the color singlet state of the QQ and F^^"^ to the anti-triplet state of the 
diquark. Corresponding relations hold for U{T) and S{T). Hence F^{T) as defined above indeed specifies 
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the free energy associated to a single static (triplet or anti-triplet) color charge in a medium of temperature 
T; Uz{T) and S^lT) denote the corresponding internal energy and entropy. The relation 

U3{T) = F3{T)+TS3{T) (10) 

can then be interpreted as the sum of the energy shift F3 (T) due to the interaction of the static Q with the 
constituents of the medium, and the energy shift T S3 (T) due to the interaction between the constituents 
of the cloud, with 5*3 (T) specifying their change in number and T their energy. It is clear that U-i{T) is 
not due to any interaction between the static Q and its partner Q, and it will thus not affect the binding 
of the heavy quark pair. 

3 The Short Distance Limit 

For the QQ system, we now have to distinguish between the attractive color singlet and the repulsive 
color octet state. In the case of the singlet, when r become sufficiently small compared to the average 
separation of charges in the medium, the QQ pair constitutes for the medium a color-neutral entitity and 
is not "seen" by its constituents. Hence for r ^ T~^, the differences Fj^^{r,T) and uj^)i{r,T) measured 
in lattice studies are simply due to the perturbative interaction of the two static charges, with no medium 
effects. In the short distance regime we thus have the temperature-independent form 

where a{r) is the r-dependent running coupling and 4/3 is the SU{i) color Casimir coefficient for 3x3^1. 

For the attractive anti-triplet case of the diquark state, the perturbative form of the direct QQ interaction 
gives 1/2 that of the singlet QQ, as determined by color SU{'i) Casimir coefficient 3x3—^3. Now, 
however, the small overall system is still colored and hence is seen by the medium as a single point-like 
color charge. It will thus again lead to a polarization cloud, which for a triplet or anti-triplet just is 
F3(r). As a result, we have in the short distance regime the free energy [6] 

F(^(r,T)^-^^+F3(T). (12) 

Comparing eqs. ((TT|l and (fT2l) . we expect that 

i^«(r,T)^2[F®(r,T)-F3(r)] (13) 

should be satisfied in the small distance limit. Using eq. ^ we see that it also holds at large distance, 
and as seen in Fig. [31 it does so for T > Tc even at intermediate separation distances. 




Figure 3: Free energy of a singlet QQ and of an anti-triplet QQ [6], see eq. ([T3| [6 
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We thus find a general pattern of behavior for the attractive case (colorless singlet QQ and colored anti- 
triplet QQ). The only caveat is that while the singlet QQ becomes simply Coulomb- like at short distance, 
for the attractive QQ antitriplet state a polarization cloud remains even in the short distance limit. Once 
this is taken into account, the remaining interaction appears to be quite insensitive to whether we consider 
a QQ or a QQ system. 



4 The Intermediate Separation Regime 

To study the behavior at finite r, we consider the effective couplings ap and ajj for the singlet QQ state, 
with 

„H.T,=!.»r2£ipi)v (14) 



A \ dr ^ 

and a corresponding relation for [/(r, T) in place of F(r, T). For sufficiently small r, there are no medium 
effects and hence ajj = ap. Since S{r,T) becomes constant for large r, they also become equal in that 
limit. The behavior of aF{r,T) for the singlet QQ system, as found in lattice QCD, is illustrated in 
Fig. [3^ for a range of temperatures T > Tc [8]. It is seen that aF{r,T) follows the vacuum form ^ 
up to the point where Debye screening begins to set in and eventually makes it vanish. In the range 
of r values of interest here, this behavior disagrees strongly with a Yukawa-like Debye-screened form 
Fi{r,T) ~ a{T)r~^ ex-p{—r/rD{T)}, with rD{T) as color-screening radius; the actual form has been 
studied in detail in [9]. As a result, models based on a Yukawa form [10] cannot correctly reproduce the 
free energy data in the crucial range of r and T; either a, or m/j, or both must depend on r as well 
as on T. In Fig. we also show the corresponding behavior of au{r^T). As argued above, the two 
couplings agree at fixed T in the small as well as in the large distance limits. In the intermediate r-range, 
the internal energy coupling considerably exceeds the vacuum Cornell form as well as that obtained from 
the free energy. The reason for this excess is that the internal energy difference contains not only the 
QQ interaction, but also the effect of the interaction of the charges in the polarization clouds. With 
increasing temperature, the correlation length and hence the size of the polarization clouds decreases, so 
that the difference between the two couplings also decreases. 
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Figure 4: Effective couplings from (a) free and (b) internal energy as function of r 



The difference between the values of U{r,T) and F{r,T) obtained in lattice studies thus becomes clear. 
The internal energy measures the binding potential between the two heavy quarks, that between each 
quark and the constituents of its polarization cloud, and that between the constituents of the overlapping 
polarization clouds. The free energy, on the other hand, measures only the potential between the "bare" 
heavy quarks, modified by the screening effects due to the polarization clouds. In the large distance 
limit, the heavy quark potential vanishes, so that F approaches with increasing temperature the entropy 
change —TS{T); since the entropy change becomes very weakly (logarithmic) temperature-dependent for 
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large T, this change grows almost linearly with T. For r — s- oo, U{r,T) approaches the twice interaction 
energy between a heavy quark and its cloud constituents; if the correlation length and hence the cloud 
size vanish in this limit, U{oo,T) is also expected to vanish. In the short distance limit, the QQ screens 
itself, so that both F and U measure only the direct QQ interaction. 

Stated in other words, the attractive interaction between two heavy quarks inside a deconfined medium 
has two distinct components: the direct Coulombic QQ interaction, and the non-Abelian interaction 
between the constituents of the polarization clouds which the static quarks acquire in the medium. At 
large distances, this "dressing" simply reduces the effective charge through screening. But at interme- 
diate distances, when the dressings overlap, they provide a strong additional attractive interaction. To 
corroborate this quantitatively, one would need to show that the color-averaged gluon-gluon interaction 
is attractive. For static quark sources of different quantum number, including adjoint octet quarks, such 
an attraction has been shown [7]. For gluons, it is not evident how similar studies could be carried out. 
The general effect, however, is conceptually similar to the energy loss of a fast parton passing through a 
deconfined medium: it loses energy through direct (bremsstrahlung) interactions with the medium [11], 
but much more through the indirect interactions of its gluon cloud with the medium [12]. 

It should be emphasized that the non-perturbative interaction observed in the intermediate r regime 
occurs (apart from the differences in the short distance limit) for a colored anti-triplet QQ state in just 
the same functional form as it does for the colorless QQ singlet. 



5 Potential Model Studies 

These considerations can throw some light on the problem of extracting from lattice results the temperature- 
dependent QQ binding potential V{r,T), to be used in a two-body Schrodinger equation, 

2me - — V2 + V{r, T)\ $,(r) = M,$,(r), (15) 

in order to study charmonium dissociation. Here nic denotes the charm quark mass, Mi that of char- 
monium state i. First studies had used the color-averaged free energy for V{r,T) [13], but subsequently 
lattice results for the color-singlet free energy became available, leading to somewhat stronger binding. 
Eventually it was argued that the correct potential is given by the internal energy U (r, T) . The proposals 
of the last years now cover the whole range, with the general form xU{r,T) -|- (1 — x)F{r,T), where 
< a; < 1 [13-17]. As is evident from what was said here, the effective binding increases with x. The 
internal energy U{r,T), as the expectation value of the difference of the Hamiltonians with and without 
the pair, includes the indirect "cloud-cloud" binding in addition to the direct QQ interaction and hence 
provides a stronger binding. We believe that the indirect binding cannot be neglected, so that the correct 
form of the potential should indeed be U{r,T). 

To clarify the situation further, let us consider the semi-classical limit of eq. ([T5|) applied to the case of 
J/ip dissociation, 

2 

2m, + ^ + U{r, T) = 2m, + ^ + U{r, T) = Mj/^(r, T), (16) 
TOc rricr 

where the minimum of the energy 

E{r,T) = ^ + U{r,T) (17) 

as function of r determines the ground state mass Mj/^ of the J/ip. The constant c arises from the 
uncertainty relation p^r^ ~ c and depends on the form of the binding potential. We determine it by 
requiring the correct J/^p mass at T = 0. From 

dE(r, T) , , 

= (18) 

we obtain 



^1 (19) 
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which leads to 

Mj^^,{r,T = 0) = 2m, + |ar - ^ = 3.1 GeV. (20) 

With nic — 1-3 GeV, a — 0.2 GeV^, and a = tt/12, this resuhs in c — 1.56 and for the radius of the 
J/ijj the value r(T = 0) ~ 0.42 fm, which agrees very well with that obtained in the corresponding 
quantum-mechanical study [18]. The minimization requirement (jl8|) for finite T leads to 

or 3 

where au{r,T) is the effective coupling determined above (see Fig. |4|). In Fig. [5l we solve this equation 
graphically. It is seen that up to some temperature value Tdis — 1.5 Tc, au{r,T) attains a peak large 
enough to intersect the kinetic term c/rricr, so that there is a minimum in the energy and hence a bound 
state. For T > Tdis, this is no longer the case, the energy decreases monotonically with r and there is 
no more bound state. The radius of the surviving J /ip is seen to vary very little with temperature; it 
remains in the range 0.3 — 0.45 fm up to Tdiss, where the bound state disappears. 




Figure 5: Semi-classical picture of J /ip binding and dissociation ranges 



The potential U{r,T) used in the Schrodinger equation has the correct T-independent form for small r, 
and enhanced attraction in the intermediate r range up to Tdiss- At large r, it leads to the r- independent 
value 2Uz{T), which, as we saw above, has little effect on the strength of the binding. Its value does, 
however, specify the binding energy as the gap between bound and free heavy quarks. With 

AS(T) = 2K + f/3(T)] - Afj/v,(r) = 2C/3(T) - ^ - f/(ro,r) (22) 

as binding energy, with rg determined by the minimization condition (|2ip. we obtain at the dissociation 
point a value of AE{Tdiss) ^ 0.2 GcV. 

For T < Tdis, the kinetic curve in Fig. [5] crosses the potential curve at two r values: the first crossing 
specifies the minimum of the energy and hence the J/ip bound state radius, while the difference between 
the second and the first provides a measure of the thickness Ar of the potential wall. As Ar — > 0, 
quantum effects (tunneling) will become more and more likely and thus provide a quantum-mechanical 
possibility of the dissociation even below Tdis, as well as a finite binding probability even above Tdis- 



6 Charmonium Flow 

Our arguments suggest that charmonium binding in a hot quark-gluon plasma is to a large extent due to a 
binding between the clouds of the medium surrounding the heavy quarks. If this medium is experiencing 
an overall flow, the motion of the clouds will be transmitted to the QQ core and cause a drag leading 
to charmonium flow. Similarly, the isolated heavy quarks (for r oo) will experience the drag of their 
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polarization clouds, and at hadronization lead to flow of open charm mesons. An observation of J/t/i flow 
in heavy ion experiments should therefore not be interpreted as evidence for primary J /ip dissociation, 
followed by regeneration due to binding of charm constituents from different collisions. The cloud-cloud 
binding discussed here can transfer any medium motion also to a primary QQ pair. 

7 Outlook 

Understanding the relation between lattice and potential studies of quarkonium binding is obviously an 
essential step in solving the in-medium behavior of quarkonia. What we have presented here is certainly 
not a final answer. It is only meant to recall 

• that the heavy quark interaction in a quark-gluon plasma is not a simple two-component problem, 
but that on the other hand, 

• this does not rule out a description in terms of a suitably formulated potential theory. 
Recent analytical work [19-22] may provide a key to further developments in this direction. 
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